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ABSTRACT 


A computer simulation of the Brayton cycle is used to develop a 
design procedure with respect to minimizing volume and weight for the 
counterflow plate-fin recuperator. Based on the Orenda OT-4 600 horse- 
power gas turbine, recuperator performance and dimensional character- 
istics are presented for an idealized equilateral plate-fin recuperator 
core matrix. The effects of plate spacings varying between 0.1 inches 
and 0.5 inches on recuperator performance characteristics are presented 
over a wide range of core frontal areas. Specific trends toward minimum 
volume and minimum weight plate-fin recuperator cores are discussed in 


detail. 
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plate thickness 
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heat capacity rate 
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mass flow rate per unit area 
acceleration of gravity 
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thermal conductivity 
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2 
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Lt 

ft 

ft 
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eve 

lbm/(hr sq ft) 
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BTU/(hr sq ft deg F) 

hp 
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BTU/(hr sq ft deg F)/ft 
ft 

lbm/ft 

lbf/sq ft 
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BTU/hr 


BTU/1bm 


English Letter Symbols (continued) 


Th 


hydraulic radius 

fin spacing 

specific entropy 

specific fuel consumption 
temperature 


overall heat transfer 
coefficient 


velocity 

volume 

solid volume 
specific volume 
available shaft work 
compressor work 

mass 


fin thickness parameter 


Greek Letter Symbols 


of 


© 


on 


ratio of total transfer area 
of one side of exchanger to 
total exchanger volume 
compactness 

fin thickness 

regenerator effectiveness 
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temperature difference 

fluid viscosity 


density 


ratio of free flow area to 
frontal area 
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deg 
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lbf hr/sq ft 
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combustor efficiency 
compressor efficiency 
fin effectiveness 
thermal efficiency 
turbine efficiency 
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PP 
Pe 
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St 
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relative pressure 

relative pressure drop 
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Stanton number 
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I. INTRODUCTION 


The heat exchanger, a device which effects the transfer of heat from 
one fluid to another, is found whenever it is necessary to transfer 
energy without mixing the fluids involved. It is felt by many that our 
future as a power rich nation may depend partially on the heat exchanger 
engineer, as the most promising sources of power often present the most 
difficult applications. One such power source, the gas turbine, has a 
number of inherent advantages over the common gas and diesel engines; 
among them are compactness, low weight, ease of maintenance, and multi- 
fuel capabilities. 

With the addition of a heat exchanger of reasonable size, the ther- 
mal efficiency of a gas turbine unit may be increased to about the same 
as, or slightly higher, than the gas engine [7].* The added weight of 
the exchanger, however, could well negate one of the primary advantages 
of the gas turbine. Furthermore, the relative rise in production cost 
due to the addition of a heat exchanger ought to be considered in any 
comparison, so as to keep the gas turbine on a specific cost level with 
the diesel engine. 

Heat exchangers, however, are inherently bulky, heavy, and expen- 
sive. To overcome this problem, the concept of the compact heat ex- 
changer has come into being. Not necessarily small, such an exchanger 
incorporates a heat transfer surface of high area density (leading to a 


high effectiveness for given weight and volume limitations) and, generally, 


* 
Numbers in brackets designate references listed in Bibliography. 
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may be manufactured at a much lower cost per square foot of heat trans- 
fer area than the not so compact heat exchanger [9]. 

Although standard lines of compact heat exchangers have been adopted 
by many manufacturers, such exchangers may not be purchased as off-the- 
shelf items, as a proper design is needed for nearly every requirement. 
Such a design procedure usually considers the factors of heat transfer 
requirements, pressure drop characteristics, size, and cost, with an 
emphasis on the balance between the gain due to heat recovery and the 
loss due to pressure drop contributed by the addition of the exchanger. 

There exist, however, numerous applications for which this much 
used design philosophy is not sufficient. Often the designer is faced 
with stringent volume or weight limitations. The technique of deter- 
mining an optimized counterflow heat exchanger with respect to volumetric 
and weight considerations is developed in this study. 

Fundamental relationships, assumptions, and design points upon 
which the investigation is based are presented in Sections II through 
IV; the computer program used to simulate the regenerated gas turbine 


cycle is discussed in Appendix B. 
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II, THE IDEAL REGENERATED BRAYTON CYCLE 


The gas turbine plant, to a large extent, must have all its units 
designed with reference to the cycle on which it is to operate. The 
designer, therefore, must have a knowledge of the particular cycle to 
be considered. The most prevalent type of gas turbine in use today, 
the simple open cycle, has achieved great success due to its versatility 
[6]. With the addition of regeneration the efficiency of this cycle may 
be greatly increased. 

The Brayton cycle, the basic thermodynamic cycle of the gas turbine, 
(refer to Figure 1) is characterized by taking air, as the working fluid, 
from the atmosphere into a compressor where part of the turbine power is 
used to raise the pressure before the air enters the cold side of a re- 
cuperator. The recuperator reduces the heat input to the cycle by trans- 
ferring heat from the hot turbine exhaust gases to the cooler air. A 
further temperature increase occurs as fuel is burned directly in the 
air flow (forming gaseous products) through the combustor. The gas then 
enters a turbine, where it expands to a lower pressure while releasing 
energy available for work. Upon leaving the turbine stage, the gas 
enters the recuperator hot side and transfers much of its otherwise waste 
heat to the incoming air, before being exhausted into the atmosphere. 

The particular cycle used in this heat exchanger investigation 
is based on the Orenda OT-4 gas turbine, a 600 horsepower unit engine- 
ered for the U. S. Navy [10]. The design specifications for the Orenda 


unit are: 


15 


le 


en 


net power output 


overall thermal efficiency 


600 hp 


34% 


Pl compressor inlet pressure 2120 1lbf/sq ft 
Tl compressor inlet temperature 560 deg R 

T4 turbine inlet temperature 2200 deg R 

P6 exchanger exhaust pressure 2120 lbf/sq ft 
P2/P1 compressor pressure ratio 4.0 

P3 - P4 combustor relative pressure drop 0.02 

a compressor efficiency 86% 

Cae combustor efficiency 98% 

ee turbine efficiency 86% 

Tey fuel inlet temperature 530 deg R 

qy fuel lower heating value 19000 BTU/1bm 


Other parameters which are usually specified for the cycle, such as the 
mass flow rate, regenerator effectiveness, and exchanger pressure drops, 
are not presented; they are, in fact, considered as variables in the 
optimization procedure to be discussed. 

A more detailed discussion of the Brayton cycle is outlined in 


Appendix A. 
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Temperature,T 
Enthalpy, j 





Entropy,s 


Recuperator 





FIGURE 1 


THE 


IDEAL REGENERATED BRAYTON CYCLE 
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III, RECUPERATOR SURFACE AND CORE GEOMETRY 


For all applications in which the turbine inlet temperature is 
limited and a low to moderate pressure ratio is to be used, the most 
efficient method of improving the economy of the gas turbine involves 
the use of a heat exchanger, in which hot turbine exhaust gases give up 
heat to the relatively cool air from the compressor. Two basic heat ex- 
changer models, the recuperator and regenerator, are in common use. The 
term recuperator (which is considered in this study) is usually applied 
to a continuous flow method of heat exchange, while the term regenerator 


specifies a periodic flow device, usually involving a moving matrix. 


A. COMPACT HEAT EXCHANGERS 

Considerations involving heat transfer characteristics of the gas 
turbine recuperative unit have led to the development of the compact 
heat exchanger. The term compact does not refer directly to a small 
bulk or weight, but to the heat transfer surface density of the ex- 


changer. The most useful measure of this density is the compactness, 


© » where 


S = total transfer area of one side of exchanger 


volume between the plates of the same side 
Somewhat arbitrarily, the working limits of compactness for a simple 
counterflow recuperator may be specified as 1000 sq ft/cu ft as the 
maximum working limit, and 200 sq ft/cu ft as the minimum limit before 
leaving the realm of compactness [9]. 
An understanding of the geometric parameters, presented in Figure 
2, and subsequent relations most often used is essential to any compact 


heat exchanger design procedure. 
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Er 
Se 
re 
b @------ 
Cool 
Air Hot 
Gases 
FIGURE 2 
COMPACT COUNTERFLOW RECUPERATOR CORE 
a plate thickness | Ee 
b plate spacing BG 
c recuperator width ft 
Th flow passage hydraulic radius oe iile 
Co) compactness sq ft/cu ft 
A heat transfer area of one side sq ft 
A. free flow area of one side sq ft 
Ar. frontal area sq ft 
L flow length et 
V volume Cl iets 
us solid core volume cu ft 
A ratio of total heat transfer 
area of one side to the total 
exchanger volume sq-ft/cew ft 
ome ratio of free flow area of one 
side to frontal area dimensionless 


The following equations give the relations between surface and 
core factors for one side of the counterflow plate-fin recuperator. 
Subscript 1] refers to any one side, and 2 refers to the other side. 


Factors without a subscript are common to both sides. 


—~ Ar —~ _Ai 
ni 7 Ware 7 Sy Lee 30K 


ea. =| ant 322 


d, = Ai a, ow oe S. 343 
V Qatar. 
= = Ac _ Airfwur — aii a Ramee. Ea 


Eine 7 ae bit Dea+ 20. 


Ae. = on Ax. par! ——— 325 


Ay — OE 3.6 


Vs =U] (Ay) -(Ag + BA "< 


Among the factors common to both sides is A the total recuperator 


fie” 
frontal area. The preceding relations are sometimes written on the 


basis that, for the counterflow recuperator, the hot and cold side 
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frontal areas are defined separately such that their sum equals Ay. 
This investigation requires no such separation of terms and is based 


on the interpretation of frontal area as presented in Figure 2. 


Bi PLATE - FIN SURFACES 
1. General Theory 

Although other matrix forms are used, the plate-fin recuperator 
matrix is the most common [6]. Not only does this type of structure lend 
itself to light weight, extremely compact designs, but the plate spacing 
on each side may be easily optimized independently, allowing for design 
flexibility. 

An understanding of the specific geometric parameters of the 
plate-fin matrix (Figure 3) and their relations is necessary to the 


design procedure. 





FIGURE 3 


PLATE + FIN GEOMETRY 


tat) 


a plate thickness 


plate spacing 


§ fin thickness 

S) fin spacing 

A heat transfer area 

A. free flow area 

PY wetted perimeter 

L recuperator length 

W fin thickness parameter 


Ton © 


Ac 


= S(b-w) 


ey 2 sy = 2(ze=,)+ i Siete + 1) 


a= ele 


ft 


ig 


re 


ft 


Sq 


sq 


ft 


ft 


Le 


Et 
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3 a2 


3.10 
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From these relationships, the values for two very important parameters 


may be found. 


Y= = = (o-w) [2 (Gore . yy 


oz 


jek 


@®=A{(4bS)L] = 2 = aleve 4) 3.13 


Since the recuperator matrix studied involves a fin angle, © , 
of sixty degrees (resulting in an equilateral fin matrix) a further 


simplification may be made. 


wane tars 3.14 
=a 3215 
> = = 3.16 


Table I represents values of these parameters for the plate spacings 
investigated. Although the spacings greater than 0.4 inches fall out 
of the previously defined compactness region (200 < = = 1000), they 
are included in the study since: (1) this region is arbitrary and (2) 
recuperator trends in weight and volume are more easily observed when 
large and small spacings are compared. 

Zo Speciilic Parameters 

In order to keep the investigation on a common basis, all cores 

analyzed are assumed to be constructed of stainless steel, with con- 
stant properties, and uniform plate and fin thicknesses. The constant 


parameters used in the analysis are: 


a plate thickness 0.0008333 ft 
(0.01 in) 

§ fin thickness 0.0004166 ft 
(0.005 in) 
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Pr 


thermal conductivity of 
stainless steel plates 
and fins 


density of stainless 
steel plates and fins 
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12 BTU/(hr sq 
ft depar)/ ft 


484 lbm/cu ft 


TABLE |! 
EQUILATERAL PLATE - FIN SPACINGS 


Plate Compactness, @ Hydraulic ~ 
Spacing,b Radius,r, 
(in) (sq ft/cu ft) (ft) 
0.1 720.00 0.0012500 
0.125 576.00 0.0015972 
G...15 480.00 0.001944 
0.175 411.45 0.0022917 
0.2 360.00 0.0026389 
0.25 288.00 0.0033333 
0.5 240.00 0.0040278 
0.35 205.71 0.0047222 
0.4 180.00 0.0054167 
0.5 144,00 0.0068056 


* 


The hydraulic radius, r is based 


h ’ 
on a fin thickness, 6, of 0.005 tn. 
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IV. HEAT TRANSFER AND FLOW-FRICTION DESIGN THEORY 


The design of a recuperator involves a consideration of both the 
heat transfer rates between the fluids and the pumping power expended to 
overcome fluid friction in the exchanger [8]. Obviously, one of the most 
desirable characteristics of a compact heat exchanger is that it yield a 


maximum heat exchange for a minimum pressure drop. 


A. HEAT TRANSFER CHARACTERISTICS 

Elementary heat transfer principles show that the cool air and hot 
gas flows should be in opposite directions; otherwise, the two streams 
would approach their mean temperature only. Subsequently, many compact 
heat exchangers involve a counterflow design, with plate-fin construc- 
tion for added efficiency [5]. 

An overall heat transfer rate equation which includes convective 
and conductive heat mechanisms responsible for the heat transfer from 
the hot gases to the cool air may be written for a differential element 


dq. =k) Qin ap) 4.1 


where U is an overall heat transfer coefficient, dA represents an 
element of heat transfer area, and dq is a differential heat flux. 
A physical understanding of the heat transfer performance of the heat 
exchanger may be attained with the use of equation 4.1 and a number of 
parameters to be defined. 

Da Fin Effectiveness 


The straight recuperator fin of constant thickness § and 
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length Q is represented in Figure 4. Such a fin extends from the 


base plate, at temperature To» through the moving fluid, at temperature 


ae , to a second plate. As there is no heat flux through the fin 


midway between the plates (due to the symmetry involved), a heat trans- 


fer analysis need only consider a fin of effective length Q /2, insul- 


ated at the tip (Figure 4). An energy balance over an element of this 


fin [2] results in the equation: 


q = km A, @o Tanh(m 2) Lyo2 
k = thermal conductivity of fin material 
ne _ = a dimensionless parameter for thin fins 
A, * §¢1) = cross section area of fin per unit width 
Sn = Ros 
h = 


convective heat transfer coefficient between the fin and 
the fluid 





FIGURE 4 


STRATGHT FIN OF UNIFORM THICKNESS 


ZI 


In the order to express the heat exchanging capacity of an extended 


surface it is useful to define the term fin effectiveness as 


fin effectiveness 


ic: 


_ actual heat transferred from fin 
rey ~ heat transferred if entire fin were at temperature Ty 


_ kmAg Se __Tanh(m 2/2) 
ae 7 hh /2) On 4.3 


(Tani ony) aGaeon/2) ia 


Ns 


Ze Surface Effectiveness 
Since the plate-fin geometry involves an array of fins, 
it is useful to define a total surface temperature effectiveness,Qo , 
which gives a measure of performance of the entire array. Defining 
"A" ag the total heat transfer area of the combined fin and plate 


surfaces, 


_ actual heat transfer from A 4.5 


Ro ~ heat transfer if A were at To 


ae Coe NTs vena Wee NIG eis ae 


A h (Ty - T a) 


Nor b> BE Ging) 


a Overall Heat Transfer Coefficient 
The overall conductance U, as defined previously, may 


now be calculated as 


im 
= ee SES eee = eee eee a7 
we | (AR/A)Ront  ~ CAr/Ad k ma 


where wy is based on a unit of cold side total surface [2]. In the 
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usual gas-to-gas recuperator, the second term in equation 4.7 (the wall 


resistance component) may be neglected relative to the fluid resistances. 


Lon overall surface effectiveness for hot and 
‘lac cold sides 
A, average wall area 
AL total heat transfer area of hot and cold sides 
Cc 
h convection heat transfer coefficients of 
r hot and cold sides 
a plate thickness 
k thermal conductivity of plate material 


on Counterflow Exchanger Heat Transfer Effectiveness 


As previously stated, the most common and most effective flow 
arrangement for a two-fluid heat exchanger is the counterflow design. 
Referring to Figure 5, an energy balance leads to a second equation for 


the amount of heat transferred between the two fluids. 

ech = Gamo a=) coven) cacdanc 4.8 
Integration over the length of the recuperator leads to the equation 
9, = Cree), Sar = ieee) = Crm Cpe CVeaeat -TVean) 4.9 


It is useful to define the term "capacity rate'’ as the product of the 
specific heat of one fluid and the mass flow rate of that flied CL 


referring to the hot gas, and C. referring to the cooler air. 


o = Cry CT Rin = OR SEC) = Cec Gi 2 soe a Tein) 4.10 
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= == a. SSeS = = =o —_—=_ lc hl — ee ee —_——— 


h,in 


c,out 





flow length 


FIGURE 5 


TEMPERATURE VARIATIONS IN COUNTERFLOW HEAT EXCHANGER 


The heat transfer performance of the recuperator may now be deter- 
mined by the heat exchanger effectiveness, € , defined as the ratio of 
the actual heat transfer to the thermodynamic maximum heat transfer 
rate. The maximum heat transfer would be obtained if one of the fluids 


were to undergo a temperature change equal to (T, . =-T ), the 
h,in Glin : 


maximum temperature difference in the counterflow exchanger. Since the 
heat lost by the hot gas must be gained by the cooler air (neglecting 

losses), as may be seen in Equations 4.9 and 4.10, the fluid undergoing 
this maximum temperature change must have the lower of the two capacity 


rate values, C : 
min 


Snax - Stan te 7 Toi a 


The effectiveness may now be determined as 


ms actual heat transfer 
maximum heat transfer 


30 


é = Cn CORRS a algae = Ga Cale out - ipee 4.19 
Cmin en An “TR in ) Comin Gen ween) 


From the previous equations, the effectiveness may be shown to be 


independent of the fluid temperatures [2]. 





Crin UA 
a Ce Glo 
aa ———— | 4.13 
c Cmen 
1 - Smin ca - Gun) A 
ie Cm aes QO. Comat af. 
where as and C are, respectively, the smaller and the larger of 


the two capacity rate magnitudes. 


Finally, defining the Number of Transfer Units, NTU (a dimension- 


less expression for the heat transfer size of the recuperator), as the 


ratio of (U A) to C in? where A refers to the total surface area on 


which U is based, the effectiveness may be written as 





Cmin 
-Ci- Sew) LNG! 
28 i - 2 —_ 4.14 
a5 
i- So Ears 


A comparison of this equation with effectiveness relations for other 
geometries will lead to the conclusion that the counterflow arrangement 
is indeed useful for gas flow recuperators (for which C /C 


min max 
when an effectiveness in the order of 80% is to be obtained 


~ 1» 
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Be FLOW - FRICTION THEORY 

In the design of a gas-to-gas heat exchanger, the pressure drop 
characteristics assume an importance equal to that of the heat trans- 
fer characteristics [8]. The designer is generally interested in the 
overall core pressure drop, including entrance and exit effects, rather 


than the core friction pressure drop alone. 
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FIGURE 6 - 7 


RECUPERATOR CORE MODEL FOR COOL AIR FLOW PRESSURE DROP ANALYSIS 


i Core Friction and Flow Acceleration 


The core friction pressure loss term, due to viscous shear 


pua 


and pressure force, is expressed in terms of the dynamic head 2 - 





the appropriate friction factor, i » and the ratio of the length to 


hydraulic radius of one side, 2 meas » as 
r A 
h Cc 
2 
aR, = ¢ A \ (Poa 4.15 
ickicon. Ac eA 
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ATction = (2) i 4.16 


where G is based on the minimum free flow frontal area (G= Pur: DaL2: SF) 
Ac 
and Om is the average core fluid density Ben = ey 
Since density changes occur in the core as the fluid temperature 
changes, a term accounting for flow acceleration must be included in 


the pressure drop analysis. The change in momentum due to this fluid 


acceleration may be expressed as 





Ee 
a Peter ota 7 (a Q 4.17 
BO GU 
fa Seeelerat Lon = ACU ae) 4.18 
Ac. Qo 
rs O 
OPoecceltration = 7 a = . 1) 4,19 


2. Entrance and Exit Effects 
The entrance pressure drop, due to the free flow area change 
as well as to the irreversible free expansion following the abrupt in- 


let contraction may be written as [8]. 





ai oe 
Pr- Rm = (2) pu? Ciuaeaey 5 joy We 2a 4.20 
Cc? 
Pi - Pa = 2D ie a- 2 en) 4.21 
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where a constant density is assumed to exist from point "1" to point 
eile Ko? the contraction, or entrance, coefficient may be calculated 
by methods explained in Ref. [8]. 

Similarly, the exit pressure rise is due to the change in flow 
area (which is identical to the corresponding entrance pressure drop 


term) and the irreversible free expansion and momentum changes follow- 


ing the abrupt exit expansion. This pressure rise is written as 





i 2 

woo cee (za) eae Co-?-4) 4 Coen —— 4.22 
Gq 

oes) eh eg meme 18) 1S)? Co: 2 We) 4.23 


where Ky the exit coefficient, is evaluated in Ref. [8]. 


33 Overall Pressure Drop 


The combination of all pressure drop effects results in the 


equation 


Ge . 2VPr Pa er A \ ) 
eee oes = 2Q Qu (leo eke) +t Ba (Be axe + Om ia » Bite) 


4.24 


As it has been assumed that Spat fos and Pe Pz the relation for 


the overall exchanger pressure drop in terms of specific volume, where 


Wy 4a V ' 
Ve = ee may be written as 





Cay V A Ven V2 s 
P,-P2 s 5g (ke a bet -20ui <4) +) evi Si -etMe) 


4.25 
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2. PARAMETRIC APPROXIMATIONS 

In the preceding heat transfer and flow friction analysis, a number 
of important engineering parameters and fluid properties are used in 
fundamental relationships. Usually presented only in tabular or graphi- 
cal form, these parameters may not be directly obtained in a computer 
analysis. To overcome this problem, functional relationships approxi- 
mating the tabulated values must be made a part of the computer program. 

For the program [1], such relationships were formulated by making a 
"least squares fit'' of the curves representing tabulated data of the 
required parameters. The resulting functional relationships, of the 


form 
Pox = Ay + A2X + Asx? dt sec+ 4.26 


where: P(x) = desired parameter 


x = variable on which the 
parameter is based 


A; = polynomial coefficients 
resulting from least 


squares fit 


i-l = order of resulting 
polynomial 


are used extensively in the computer analysis. Table II specifies the 
parameters for which such polynomial approximations are used, the 
variable on which each is based, and the order of the polynomial. 

Table II indicates that two parameters for which functional relation- 
ships must be obtained are dependent on the Reynolds number, 04 (40 )/30 
The Fanning friction factor, n ,» is used in equation 4.15 and subsequent 
relations to obtain the pressure drop due to core friction through the 


recuperator. The Colburn-j factor, j, has been obtained experimentally 
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for a number of recuperator geometries [8] and is used to obtain the 
convection heat transfer coefficient, h, as follows: 


This heat transfer parameter is defined as 


‘ur Seabee 4.27 


Since j and Pr are determined from functional relationships, the 


Stanton number may be calculated as 


st = 4 Pre! 4.28 


Thus, the convection heat transfer coefficient, used in equation 4.5 


and subsequent relations, is obtained as 
pe he Tee u 4.29 


Experimental investigations of t and j provide curves similar 
to those presented in Figure 7. As indicated, the laminar flow regime 
(Re < 1000) produces linear relationships for the ideal equilateral 


triangular passage such that 


rh 
it 


13.33/Re 4.30 


2.64/Re “ime 


Need © 
li 


The turbulent region (Re 7 5000) data may be closely approximated by 
the linear relationship resulting from circular tube flow as 


0.046 Reon 432 


rh 
il 


0.023 Re 9"* 4.33 
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FIGURE 7 


EXPERIMENTAL CURVES FOR { AND j 
The transition region (1000 £ Re £ 5000) is also based on experi- 
mental data and must be approximated by a polynomial equation for 
computer analysis. The polynomial coefficients for the resulting 
equations approximating this transition region, as well as those for 


the other parameters in Table II, are presented in Appendix B. 


37 


TABLE 1! 


FUNCTIONAL RELATIONSHIPS 


Parameter 
Fanning friction 
factor 
Colburn-j factor 
Prandtl! number 
viscosity 
specific heat 
enthalpy 
relative pressure 
temperature 


temperature 


Function Of 


Reynolds number 


Reynolds number 
temperature 
temperature 
temperature 
temperature 
temperature 
relative pressure 


enthalpy 
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Order Of 
Polynomial 


10 


10 


V. RESULTS 


De, OBJECTIVES 

For many years engineers have designed heat exchangers by first 
assuming core dimensions and then calculating mass flow rates, heat 
transfer rates, and pressure drops which would occur if fluids were 
flowing through the assumed core. Extensive investigations of various 
heat transfer surfaces offer the designer enough latitude in surface 
geometry to successfully meet the heat transfer and pumping power re- 
quirements for a specific application. 

Texts frequently refer to the possibility of optimizing the inter- 
nal core dimensions of a heat exchanger once a particular surface is 
specified, yet very little research has been published in this area. 
Such a technique forms the basis for this study. 

The computer program [1] used in the study simulates the open gas 
turbine regenerated cycle under constant load conditions. All adjust- 
ments in internal exchanger geometry are reflected directly in terms of 
the cycle thermal efficiency and the overall core size necessary to 
obtain that efficiency. As is outlined in Appendix C, any computer 
simulation of the gas turbine cycle requires set tolerance limits on 
various cycle design points which must be met. Early investigations of 
general recuperator trends in the study involved a tolerance limit of 
3.3% for the desired thermal efficiency of 0.34 and a limit of 0.2% for 
the desired horsepower output of 600. More detailed analysis (as pre- 
sented in Figures 15 through 22 and Tables V through VII) required 
tighter tolerances of 1.45% for the thermal efficiency and 0.09% for 
the horsepower. The use of two tolerance levels slightly affects out- 


put data for the same core, as may be noticed. 
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B. GENERAL RECUPERATOR CHARACTERISTICS 

Trends in external dimensions of the plate-fin recuperator due to 
changes in internal spacings may be noted by investigating a number of 
different cores. A sufficient number of cores may be obtained by pair- 
ing the plate spacings presented in Table I. For convenience a parti- 
cular core will be referred to in terms of the ordered pair (bb, ) where 


bo represents the cold side spacing and b, represents the hot side spac- 


h 
ing of the core. 

Due to addition of fuel in the combustor, the mass flow rate through 
the recuperator hot side is larger than the cold side flow rate by a 
factor of the fuel to air ratio, Cela The hot side mass flow rate must 
be accommodated, then, by an increase in fluid density, flow area, or 
fluid velocity (m= e Av). Since fluid density decreases as tempera- 
ture increases the engineer is left with only two alternatives. 

A velocity increase will not only accommodate the increased mass 
flow rate, but it will also increase the heat transfer rate per unit of 
surface area by a factor something less than the first power of the 
velocity. The friction power expenditure, however, will increase by a 
factor between the square and the cube of the velocity [8]. It is this 
effect which serves as the basis for the design restriction that the hot 
side plate spacing never be less than the cold side spacing. 

General recuperator volume trends may be noted in Figure 8 for five 
specific cores. Each point on every curve represents the external core 
dimensions which will satisfy the specified design points. As might be 
expected, from the definition of the compact heat exchanger, larger cores 
are required as one moves from the more compact to the larger spacings; 


for any one frontal area, more length, and thus more heat transfer area, 
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is required as spacings increase. Excluding the extremely compact core, 
(0.1, 0.1), points of minimum volume for each core fall below a frontal 
area of five square feet. Furthermore, every core requires a frontal 
area of at least two square feet before the specified cycle design points 
are met. 

A more detailed view is presented in Table III, in which the minimum 
volume of the (0.1,0.1) core is seen to fall within a frontal area range 
between 30 and 40 square feet. Definitely, the minimum volume of such a 
core does not present the optimum package. Similar characteristics were 
noted only for two other extremely compact cores, (0.1, 0.125) and (0.125, 
0.125). 

An increase in length over a region of increasing frontal area is 
noticeable in all cores except those having a cold side spacing of 0.1 
inches. Occurring at a frontal area greater than that for the minimum 
volume core, this effect becomes more apparent as core spacings increase. 
Table III indicates that the mass flow rates for these cores decrease as 
frontal area increases. Not shown is the fact that the recuperator hot 
and cold side Reynolds numbers also decrease with the mass flow rate. 

The first knee of each curve (that point at which L begins to increase 
with A.) Occurs when the two Reynolds numbers approach a magnitude of 
2000. Reference to Figure 7 indicates a decrease in the Colburn-j fac- 
tor with decreasing Re at this point. Evidently, the resulting decrease 
in the convection heat transfer coefficient, h (which is proportional to 
j), is accompanied by an increase in heat transfer area (as L increases 
with A.) in order that the desired cycle thermal efficiency be main- 
tained. This effect continues as j decreases until the second knee at 


which the two Reynolds numbers have approached a magnitude near 1000. 
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At this point, since j begins to increase with increasing frontal area, 
less heat transfer area is required, and L begins to decrease. 

The effect is not noticeable for those cores of low plate spacings, 
as both the hot and cold side Reynolds numbers of such cores fall in the 
laminar region for all frontal areas. Reynolds numbers for cores of low 
cold side spacing and large hot side spacing (i.e., (0.1, 0.5)) however, 
tend to have a maximum Reynolds number near 2000 at low frontal areas, 
followed by a rapid drop as frontal area increases. The slight length 
increase required by such cores is not noticeable due to the magnitude 
of the length increments used at this point of the investigation. 

Relative pressure drops on both sides of all cores decrease to zero 
as frontal areas increase. In the region of minimum volume (ApS 5 sq 
ft), the cold side relative pressure drop appears to be negligible for 
all cores except that which combines a small cold side spacing and a 
large hot side spacing. The hot side relative pressure drop, on the 
other hand, appears to be high only for those cores in which the two 


plate spacings are of equal magnitude. 


or MINIMUM VOLUME AND WEIGHT CORES 

Figure 9 indicates dimensional requirements for various cores in the 
minimum volume region. As is to be expected, required lengths, and there- 
fore volumes, increase as core spacings increase. The effect of forming 
cores of unequal spacings as compared to equally spaced cores is presented 
in Figure 10. 

To this point a trend has developed in which smaller spacings re- 
quire smaller cores. Based on this trend, one might incorrectly reason 
that the minimum volume for a given hot side spacing may be obtained 
simply by decreasing the cold side spacing to a minimum allowable magni- 


tude. In fact, for the cores of 0.25 inch hot side spacing this appears 
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to be the case as b. is reduced from 0.25 inches to 0.1 inches. How- 
ever, the trend presented by those cores for which bE is 0.5 inches as 
i: decreases should be noted. A decrease from a cold side spacing of 
0.5 inches to 0.25 inches reduces core volume; a further decrease to a 
bo of 0.1 inches produces a core of extremely large volume. This core, 
which combines the large length and frontal area requirements of the two 
range limiting spacings, results in the largest minimum volume point re- 
corded of all cores investigated. 

Similar analysis shows the effect, presented in Figures 11 through 
13, of increasing the hot side spacing for a constant bo: For all cores 


presented, the b, increase results in increased frontal area and length 


h 
requirements; thus, the minimum volume for each series of constant cold 
side spacing cores increases as the hot side spacing is increased. 
Recuperator dimensional and performance characteristics at minimum 
volume points are presented in Table IV; each set of cores listed con- 


sists of a constant bo with increasing b The optimum core with re- 


h’ 
spect to volume in each set occurs at or near that core configuration 
for which by is twice bee The optimum core with respect to weight for 


each set consists of the same, or a slightly higher, b, as the optimum 


h 
volume core. Such a trend might lead one to believe a rule of thumb 
often used by manufacturers that the core hot side spacing be twice the 
cold side spacing. Analysis will show this to be an inadequate design 
procedure. Further investigation in this area, however, should be pre- 
ceded by mention of the performance characteristics of the minimum core 
points as tabulated in Table IV. 

For any one set of cores, the cold side relative pressure drop in- 


creases with increased b while the hot side relative pressure drop 


h’” 
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decreases. This effect is presented graphically for one set of cores in 
Figure 14. The large magnitude of the hot side relative pressure drop 
for cores of bE a ba is to be expected since the hot side mass flow rate 
must be accommodated by an increase in either flow area or velocity over 
that of the cold side flow. Since there is little or no area increase 
(bd, x b> velocities on the hot side of such cores must be higher, 
yielding relatively large pressure drops. As mass flow rates generally 
fall and flow areas increase as bE is increased for each set of cores, 
the hot side relative pressure drop falls. Evidently, the cold side 
relative pressure drop is affected more by the general length increase, 
as bE increases for each set, than by the decreasing trend in the mass 
flow rate. Thus, the cold side relative pressure drop rises with bE in 


each set of cores. 


D. OPTIMUM CORES 
Concentrating once again on dimensional aspects, Figures 15 through 


19 represent the effect of increasing b, for a constant oe in the region 


h 
of optimum core spacings with respect to volume and weight. It is ap- 
parent from these plots that an optimum design procedure is not that for 
which bE = 2b: Tables V and VI present the characteristics of cores 
optimized with respect to volume and weight; Table VII, which consists 
of performance data for cores of bE = 2D» is presented for comparison. 
A comparison of Tables V and VI indicates that the minimum volume core 
for any one Dy involves a higher mass flow rate and less length than the 
corresponding minimum weight core. This results in a lower cold side 


and higher hot side relative pressure drop through the minimum volume 


core as compared to the minimum weight core for any set bee 
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Further investigation of the Tables indicates that the minimum weight 
core for any be requires a larger DL than that for the minimum volume 
core. This effect (noticed earlier in Table IV) is to be expected, as 
the minimum weight core represents a compromise between a low core volume 
and a low overall core density (which decreases as spacings increase). 
Differences between the two sets of optimized cores and those cores for 
which b = 2b . are illustrated in Figure 20. 

Allowing a tolerance limit of 0.004 inches either side of the optimum 
hot side spacing for any given cold side spacing as presented in Tables 
V and VI, smooth optimum spacing curves may be drawn, as shown. In the 
range of cold side spacings between 0.1 and 0.225 inches a linear re- 
lationship closely approximates each of the two optimum spacing curves, 


as 


b, = 0.86 b + 0.12 (O.1 £ b << Wes) oa 
h c @ 


for optimum volume spacings and 


b, = 1.16 b + 0.12 (Ghleeob. < 08225) — abe 


for optimum weight spacings. 

Outside this range of cold side spacings the optimum core plots 
tend to level off, thus tending to keep the cores within the limits of 
compactness of the compact heat exchanger region. 

Since the engineer is usually at liberty to choose both plate spac- 
ings, he must have an idea of the relative volumes and weights between 
the optimum cores which he may use. Figure 21 presents volume character- 
istics of optimized cores from Table V with respect to the minimum volume 
core (0.125, 0.225). As may be seen, if the engineer expects to use a 


cold side spacing equal to or greater than 0.175 inches, he must accept 
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the fact that his optimum core will have a volume at least 20% greater 

than the minimum available. Similarly, Figure 22 represents the weight 
characteristics of cores from Table VI with respect to the minimum weight 
core of (0.175, 0.332). Although the variance in weight is not as great 

as that for volume, on a percentage basis, it may be seen that the engineer 
must expect weights at least 5% greater than the minimum available if he 
wishes to use cold side spacings other than those between 0.125 and 0.235 
inches. The very large increase in relative weight as extremely compact 
cores (b. = 0.1 inches) are optimized should also be noted, 

For the design points of this study, if the engineer were free to 
optimize any non-dimensional parameter (as mass flow rate or relative pres- 
sure drops) while keeping within 10% of the minimum volume possible, refer- 
ence to Figure 21 and the use of equation 5-1 would indicate a freedom of 
ba choice between 0.1 and 0.15 inches and bE choice between 0.206 and 0.28 
inches. A similar design procedure with a limit of 5% over the minimum 
weight core (using Figure 22 and equation 5-2) would yield a b. range be- 
tween 0.13 and 0.225 inches and a DL range between 0.27 and 0.38 inches. 

A reference back to Table IV will indicate that the minimum weight and 
minimum volume cores for all cores investigated do fall within the indi- 
cated ranges of spacings. If both volume and weight of the core are 
limited within the specified 10% and 5%, the range of allowable cold side 
spacings is further reduced to between 0.13 and 0.15 inches, and the allow- 
able hot side range to 0.27 to 0.28 inches. 

Although the design limitations outlined are applicable only to the 
specific core geometry studied, the procedure may be applied to all simi- 
lar geometries. As indicated, when core and weight tolerances are speci- 


fied, the engineer does have some freedom in choosing internal core 
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dimensions based on other requirements. In all applications for which 
the allowable volume or weight of the gas turbine is at a premium, it 

seems that the engineer would be wise to follow a design procedure for 
the heat exchanger, as outlined, which optimizes the limited quantity, 
before moving on to further core optimization with respect to perform- 


ance. 
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FIGURE 15 


FOR b, = 0.1 INCHES 






























elo five tat tea tet THEE HIE bid dest dite lecdi cate sssafaiet 
= HUGHES EEE EL HE a ce ane ieaaae 
ee ae a 
Raitt Ee 
see LEE SEER re et ee ee ue ETS 
fiat eetttntiit EpbdEEZIZED| ceed eee eee aaee cee eet add avd even ected EPEC Httcdze 
Hu EEE fe eee 
ee ae eee 
ee ee 
ppSee Sues SEEeES ceece -GUEEE CUTE SETST Sooo Gouee COUGECUGGH GeEEEooee GeEECGECCE Suree cucee CoCSeCCECE STtz rece: sescecc 
Soe ee eee eee eee 

Seaut soseut veut tes tutes seosesetan spesetenercveneqencrsntatstess’ Spsus seuss cusueuuuns suescuusesessenlse 
CGR el ge eS Fen ged EEE EEE SEE EEE 
LIES ase UE HE EE tHE GEESE EE dl SEEPS ELSE 
SE gt dae 
Sra acre reece et ee ee eee ae ea ee een, eae eee EEE 
Fee Bsa? socal sadtay gtuz seaseareas estates Fret Sse Coser estes seeeran cs jfeem asian teed Weafaee ieeeoae 
yet Scccsunseesneecceeesee seoseee sescsecuuee Go legsucecacseqss sfice Seoenut te seecenesee seeceguane seeceesenecccccans 
Reap Geeen Geseeseyee Jesee sees seeeeeeues seeeue eee Ee Gene Goeee Geeee Gee. 4 Geesceeeee coeeceseun Geeceseccacanouane 
ae EFSF SSeeHUTE PESEEEEVe ence EO See DOE eee cee cent ceed nrevOn i. Stseese =feses “S55 teed cstecissstase 


Saas fe *)d Peete assess Setcaue seeuesseeseecsscuses sesceSeugesueeessuceauecsauee. -. ‘Seaneessasee 0-000 seeeesue00 S0en5uee 
‘een ad aan anes ae Geeaa seeenn ggoao suees dacagga Ga aaa Ba ces nen ‘2 CSEeghaseacce secee coo Fane ssuee.-ae5 J agaga Bagge oaB 
sage sewsesuees sueccsucecse HH Gann agee3 seseas HH suecencusesnece auegssecs anes. Sabn Clans i. 1a Het 


: ; ceeae i} HH DESH O€eRaes ae a OES BEDS CERCR FOR STATA ES PURRR Faees Ganda oe oot ‘aan GESGeSGN008 A sesee 
seis saris a SES Eas dsere eat t far need cent tasanredeointae erect ty ttarato oe eee 


ae sect jetted eet eee Hear ant} peal fse tel ei as sez, Saesraeee 
a 


an “ot 


EE aera eye ee ma ee 
Sle epessose eaaatttt et ROH neta cones danasrey “agECEE TE 





eean 
Ms a 
‘oe 
a 
| 
oe 
i 
a 
a 


| 
Log To oo 
Sul BO SS888 
BGG S54 


hig 
‘ H 
as808 
oun 

aorae 

2 

a 

: 

2 

2 

a 


a 
f sas souassasas sasaasases sess aeressesetsehsns incite ofa teeaeail 
ea a a en gece ea ee ea 
sees GESEs couee cede Cases Genes Gnees caee cases crest cecee GeGEBRSTSh, cy TS ESBSCS. 
seus cease cuss sense seuss segue Guage cueed GUsEESCESGSEGSSSEECSGSGESE. o--cReEs) JsescusescEsesEEE 
BSEEaceeae fezee Cer Statist cpr doseeeneg soeeeeees 
a Sebet eed retee ceeds cbets seers epcee setae feces cbend seeee sees 
Breer aa ee | | WY ry -* €88 8800 Sanag auRa 


aaa goege H HH SrrroHte = P| = oBEISEEeS Seeag8 Cea 
Bee sbsos casteeSes sestes 


: pisces seesteese 
ro seers Secaeeecese sageecoas seeei tess suit: dee anaseos Sosneci 


er Gsa0eas Pitas te jit 


E ti ,/ sh) = 


cr er 





SSS ered cues oct bts tues faved od ose de eet 
ae 


GEG 2088 (e888 HHH HH 


eta ihe 


seese Bal aasrfeeed cece ScSecceteeccsstee HERA 
eases aeee GenegeGsaca senecceesa eaa egaasa 


SESrcrecs Sees CEEEs C4SES C0cEs LEVEE SEESE EEEESEEEES sonea stsasenesese ae 


GO80 1098 8 GGGSE GREER NG GES OUORA 108 Of ORNS8 OBSSOSAGSR0R8 @eaaoa 
ef pf ff 
ScCSseCessacenescerssserss cossacees CoeseGesesceseseies Ceeseeseae Genet 
eo eases een ea fp sasencedeaa gabe 





Sal 
+30 
AES, 
78S); 
29 
Zee 
cael 


(jabs joes y 


61 


3.4 3.6 Bae: &.0 4.2 4k 4.6 4.8 
L (ft) 
EFFECT OF INCREASING by, ON CORE DIMENSIONS IN OPTIMUM REGION 


aa 2 


3.0 


FIGURE 16 
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APPENDIX A 


THE REGENERATED BRAYTON CYCLE 


A. THERMODYNAMICS 

Ideally, the basic open cycle consists of the four processes of 
isentropic compression, constant pressure heat addition, isentropic 
expansion, and constant pressure heat rejection [4]. The ideal regener- 
ated cycle differs only in that the heat addition process is accomplished 
through the use of a recuperator as well as the standard combustor. With 
these points in mind, the ideal cycle (refer to Figure 23) may be des- 


cribed using the equalities: 


Poa = Pasa = Oia A.l 
Poa = Poa A.2 
To, = Téa A.3 
T3, = Teo A.4 


The most obvious differences between the actual and ideal cycles, 
from the thermodynamic viewpoint, are the increases of entropy and drops 
in pressure due to actual working conditions. These differences may be 
best explained by taking the cycle in its respective stages. 

i Compression 

State point "2" lies to the right, at an increased entropy, of 
ideal point '2a'' due to the irreversibility of the compression process; 
both points, however, lie on a line of constant pressure, since the 
compressor outlet pressure is specified by working conditions. The 


resulting equality is 


rg 






Temperature, T 
Enthalpy, i 


Entropy,s 


("a" indicates [deal 
state points) 






ae. 
W 


Turbine 


FIGURE 23 THE ACTUAL AND IDEAL REGENERATED BRAYTON CYCLES 
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a Heat Addition 
Two changes are apparent when one compares cycle conditions at 
the recuperator exit point. First, the actual pressure level is below 
that of the ideal, due to the pressure drop as the air flows through 
the recuperator. 


Peo A.6 


Second, the temperature of the actual case is below that for the ideal, 


since less heat is transferred than for the reversible process. 


< 
T, T,, A./ 


In order to raise the temperature further, the air flows through a 
combustor. The inlet temperature of the turbine is specified by its 


material properties, 


eee b : A.8 

4a 4 
However, point "4" does lie at an entropy state greater than point '4a" 
due to the irreversibilities involved between the compressor and the 


turbine. Furthermore, an additional pressure drop occurs as the air 


flows through the combustor. 


Se Expansion 


An additional entropy increase results from the expansion 
process in the turbine, moving point "'5'' to the right of the ideal state 
point. As the turbine outlet pressure is not specified, and a further 


pressure drop is expected as the gases pass through the recuperator, 
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the pressure at point "5" is necessarily greater than that for the 


ideal cycle. 


jy ad A.10 


4. Heat Rejection 


As indicated, the pressure drop experienced by the hot gases 
in passing through the recuperator results in the pressure at point ‘''6" 
being equal to the ideal pressure, since this point is considered to be 


the exhaust to atmospheric conditions. 


P 


ae Pe = 14.696 psi A.1l 


Finally, since less heat is transferred from the hot gases to the in- 
coming cold air as is ideally possible, the temperature at point ''6" 


is greater than that for the ideal case. 
TT 
Te 7 TE, AwTZ 


B, CYCLE "SYSTEMS 

An understanding of the efficiencies of each system is necessary 
if one is to work with the Brayton cycle. This may best be accomplished 
by Looking at each system process individually. It must be noted that 
the assumption of steady flow and no potential energy change of posie 
tion is made for the cycle. More detailed information may be obtained 


in Refs. [3] and [5]. 


i Compressor 


R« = compressor efficiency 
e ¢ e 1 -j 
work of isentropic compression = la 2a 
Qe = : —_—_—_—— 
actual internal work iy “i, 
A.13 
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2 Recuperator 
& 


recuperator effectiveness 


c actual heat transfer 
maximum possible heat transfer 


= CA Cpa CVn in -Wjout) anil Cn Cee owt “Vern) 14 


G Map ein CT, at. = aac ) On Cpdenun CT in - ten 


A more detailed explanation of recuperator effectiveness is 


presented in Section IV. 


sip Combustor 


80; = combustor efficiency 
OMb 


= actual energy received by gases 


ome energy released during complete combustion 
; GRE a : 
Rooms = Gs+Pryo \ie -1a - Ly ip A.15 


PS io O4% 


This definition may be more easily understood when an energy bal- 
ance based on one pound of air and L¢jgq pounds of fuel per pound of 


air entering the combustor is made. 


Vile i¢ + Cle ony + {a = Giverny 14 + See ayy, Qe A.16 
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FIGURE 24 


COMBUSTOR ENERGY DIAGRAM 


of the three entering energy terms, it can readily be seen that i, (the 


enthalpy of the incoming air) and (the fuel to air ratio times 


Tefait 
the enthalpy of the fuel) represent molecular internal energy and flow 
work. The third term represents the chemical energy stored in the fuel, 
eee since complete condensation of water vapor in the air flow will 
not occur, the lower heating value of the fuel (which disregards the 
latent heat of condensation), qd is usually used in gas turbine anal- 
ysis [3]. 


The fuel to air ratio, a very important parameter in the gas tur- 


bine analysis, may now be shown to be 


AW 
ee Turbine 
M4 = turbine efficiency 
ON a actual internal work 
work of isentropic expansion 
Xe = ee ee A,18 


laa - Lae 
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Oe Additional Calculations 
Three final calculations are frequently useful in a complete anal- 
ysis. 
Overall Thermal Efficiency 


Qa = overall thermal efficiency 


a Be net work 
Un heat supplied 


_ (hi Ptie Cig- is) - Cie - i) - 
Horsepower 


Using the conversion factor 2544 BTU = one horsepower- 
hour, the equation for the horsepower output may be seen to be 


HP = Qinls/a)CQ s)G) A.20 


2544 


Specific Fuel Consumption 


fuel burned 
net work 


SFC = (Peja) Con) A.21 


HP 


SFC 


src = 24 A.22 
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APPENDIX B 


COMPUTER PROGRAM ESSENTIAL FEATURES 


A GENERAL ASPECTS 

The computer program referred to in Ref. [1] has been greatly ex- 
panded by Professor P. F. Pucci since its original use, and may now be 
used to investigate parallel and crossflow, as well as counterflow, 
recuperator characteristics. This counterflow analysis is based on 
design parameters for the Orenda OT-4 gas turbine only; the program, 
however, is applicable to any gas turbine recuperator analysis, as the 
desired design points are introduced as input data. A basic under- 
standing of the program logic may be obtained by referring to the simpli- 
fied flow chart, Figure 25; a sample program output is presented in 
Figure 28. 

The program is written on the basis that changes in recuperator core 
geometry (plate spacings and fin configurations) and dimensions (frontal 
area and length) are reflected directly in the overall cycle thermal 
efficiency, rather than in the often used design variables of heat ex- 
changer pressure drop and heat transfer characteristics. Although the 
program may be used to obtain data for a number of conditions (such as 
the computation of: (1) the maximum cycle thermal efficiency given a 
recuperator of specified volume, (2) recuperator dimensions necessary to 
obtain a specified thermal efficiency and, (3) general outputs for a 
simple open cycle, to mention only a few) only two output types were used 


extensively in the analysis. 
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The first, Output Type III, computes the required recuperator 
length for a specified frontal area and internal geometry; furthermore, 
the horsepower and cycle thermal efficiency are specified, within pre- 
determined tolerance levels. This output is extremely useful in provid- 
ing an overall picture of the desired recuperator’s behavior over a wide 
range of frontal areas. Information on the effect of changes in internal 
core dimensions, may be obtained with the use of Output Type VIII. This 
routine, similar to Output Type III, computes the required length of a 
recuperator over a range of frontal areas. In this procedure, however, 
the cold side plate spacing remains constant, as before, whereas the hot 
side plate spacing is incremented over a desired range. Flow charts for 


these routines are presented in Figures 26 and 27. 


B. PARAMETRIC APPROXIMATIONS 

As explained in Section IV, the computer program makes use of poly- 
nomial approximations to obtain a number of fundamental parameters. The 
polynomial coefficients used in this analysis will be presented, should 
a need for them arise in the future; unless otherwise indicated, the 


polynomial is of the form 


P(x) = A, + A,X + a,x’ tT” 5 65 eaenerene tee et raeierenat: Ax? Bel 
where: P(x) = desired parameter 
x = variable on which the parameter is based 
A. = polynomial coefficients 
i-l = order of polynomial 


no 


Fanning Friction Factor 


(function of Re; 


A, = -1.311518618929068 x 10 ‘te 
A, = 7.522772995845896 x 107* ee 
A, = -1.550552244236629 x ie Aa 
A, = 1.730155685156651 x lon on 
A, = -1.189041200805969 x 1? AL = 
A, = 5.329916594794135 x 10°1° 

Colburn-j Factor (function of Re; 1000 Re 
A, = 1.644630125351775 x 107! A, = 
A, = -6.796500408006909 x ee er 
A, = 1.238993655152145 x lone aves 
A, = -1.292429783097030 x lai A, = 
A. = 8.543160388163581 x om: A, 
A, = -3.740734237674230 x iia 

Prandtl Number (function of temperature) 
A, = 0.7865630968 Ae 
A, =-0.2871584338 x 1e*" Loe 
A, = 0.3666074673 x onc Aa = 
A, =-0.1373281188 x 107? — 
A. = 0.1569904076 x as eae 

Viscosity (function of temperature) 


JA = 


14.822 


1000 < Re + 5000) 


-1.590750036125343 x 107'" 
3.138452606608756 x 10°-> 
-3.933692644883426 x 107°’ 
2 .838094432805311 x 107>* 
-8 .976537263266828 x 107>° 
< 5000) 

1.100442463486636 x 107)? 
-2.152007282686803 x 1077> 
2.683333339022575 x 10°~/ 
=141930653368061576 1s 1omeu 
6 .099637443396928 x 107>° 
-0.1973361281 x 10°) 
61041204710 «07 ~’ 

-0. 2047252924 x 1077! 
sol 1sseg01ea1 x 10m 
0.7960530427 x 1077" 

B.2 


(Q) x (1.716 x i) x 3600.0 


where Q is determined by the polynomial equation 
(as a function of temperature) of coefficients: 
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A, = -0.4132588012 x 10 A, = -0.5752187450 x 10 

A, = 0.2706898585 x 107 dees) 0.2117993215 x son 
A, = -0.1318733093 x 107° A, = -0.3862860359 x 107~* 
A, = 0.7538864121 x ry jane Ay = 0.2782105119 x Tore 
A. = 0.7221455415 x ion Ay o= 0-6763288488 x ior“ 


Specific Heat (function of temperature) 


0.24913875 x 10° A, = -0.21358038 x 107-9 


> 
tl 


- 4 
= “ _ -15 

A, = -0.48196782 x 10 A, = 0.48224330 x 10 
A, = 0.68100588 x 107’ A, = 0.47985777 x 10718 


Enthalpy (function of temperature) 


A, = -0.97746770 x 10° A, = 0.30219636 x om 
A, = 0.25257199 x 10° A. = -0.94272843 x 10° "* 
A, = -0.31167200 x 1one A, = 0.11868992 x [on 
A, = -0.32615268 x 107°” 
Relative Pressure (function of temperature) 
aa = Oat eS 9eD) Bas 
where D is determined by the polynomial equation 
(as a function of temperature) of coefficients: 
A, = -0.4043080 x 10" A, = 0.21127191 x 107 
A, = 0.2759446 x 107" A, = -0.67216409 x 107"! 
A, = -0.45085958 x 10" A, = 0.12101712 x 10°" 
A, = 0.53019327 x 107/ A, = -0.94003173 x 107°? 
A. = -0.41468032 x 10° °° 
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Temperature (function of enthalpy) 
The temperature is computed in the standard polynomial form, as a 


function of x, where 


x = log(10.0 + temperature) B.4 
A, = 0.23796617 x 10° A, = -0.45916847 x 10° 
A, = 0.37807000 x 10+ A, = 0.33772184 x 107+ 
A, = 0.49211414 x 107 A, = -0.13466334 x 1077 
A, = -0.15095696 x 107 A, = 0.23704971 x Tome 
A. = 0.36226837 x 10/ 
Temperature (function of enthalpy) 

A, = 0.40000290 x hoe A, = 0.86703630 x Om 
A, = 0.41827332 x 10! A, = -0.27042766 x TO 
A, = -0.64799748 x 107? Ag = 0.48670442 x er 
A, = 0.12160439 x 10~* Ay = -0.47242455 x TO 
A, = -0.15484329 x its Aj, = 0.19149538 x ones 
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parameters 





Calculate thermodynamic properties at cycle 


points 2 & 5, & Pe le & P 


a cu ate 
ii, mallee a 


Recalculate 
recuperator pressure drops 


Is new P_ within 


tolerance level of old 


P,? 



















Yes No 
Recalculate cycle parameters 
based on new P. 
Calculate en & HP 
is HP 
decreasing? 
Yes No 
its HP within 


"Destred HP 
Not Obtalinable' tolerance level of 
desired HP? 


No 


Calculate new hth 
from present 
Darameters 

Is Tan within 
tolerance level of 
desired in ? 


FIGURE 25 PROGRAM FLOW CHART 
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Keep Agemax’ t1? bmax’ 
desired Qn, A&al 
Increments 


Ls Ly 


Calculate outputs In 
main subroutine 
is HP within 


tolerance level of 
desired HP? 






Yes 


Increment L 


Is Yun within 
tolerance level of 
desired fin ? 


Yes No 
Is Tw decreasing? 
Yes No 


‘run unobtainable! IisMu too low? 












Yes No 


Is L too large? 


Decrease L 


Yes 


i oe 


ax 
[Print out puts| 


fr f rmax 
Yes No 
(Stop ) 
FIGURE 26 OUTPUT ttt FLOW CHART 


by Obmax? Oh 


Increments & 


Output II! Inputs 


ce) 












Compute Gy & Fr, 


Print Gy, & ip 


Go to Output III 


Is b b ? 


Yes 


Increment b, 


FIGURE 27 OUTPUT VIII FLOW CHART 
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SUMeUNS VYPE THREE 


SPECTAL INPUT PARAMETEFS 


MINIMUM 
MAXYMUM 
FRONTAL 
THERMAL 


FRONTAL AREA = 
FRCATAL AREA = 
ARE INCREMENT = 
EFFICIENCY =(.340 


2-309 SQ.FT. 
20509 SQ.FT. 


0.905 SQ.FT. 


OUTFUT CESIREC - CCRRESPONCIANG LENGTHS 


HOR SePOWER 
5399292 
TRANSFER 
UNITS 
0 UG 


3 
AIR SIDE 
DELTA P/P 


x,CG9I 
rg 
aiel rash, 


FRCNTAL AREA = 2. 


HCRSEPOWER 
Boo | 
t 
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LMIN = 
LMAX =20.0J00 FT. 
LINC = 0.095 FT. 


oie 350 FT. 


MAL EFFICIENCY NOT ORTAINABL ITH FRCNTAL AREA = 2.3CG 
MAL EFFICIENCY =—-9.337 AT Saye fe GBF S.3508 i. 
TREF MAL SeSCI FIC FUEL 10 AIR MASS FLC 
SEFICTENC FUEL aCON RATIO RATE 
baye  Salee O.397 0.C108 22140. 
REGENERATOR AIR SIDE GAS SINE 
CMIN/CMAX EFFECTIVENESS REYNOLDS NC REYNOLDS A 
42966 7. 759 4109. 371 Sie 
GELS SIDE REGENERATOR RFGENEPATCR ATIC GAS/ATI 
DELTA P/P VOLUME WEIGHT FRONTAL ARES 
a ae 12.3050 399.71 1.292 
P 3 P+ P5 Nez T3 1S T6é 
Poke oat re |) 15.75 B73. 1487, 1681. 10 38. 
205 SQOsm ie REQLIRED LENGTH = 5.680C FT. 
THERMAL SPECTE TG FUER ST CeArR MASS Els 
EFFIC ] GimMCc Y FUEL “CEN RATIO RATE 
weal <a 32394 (0 eG106 222759. 
REGENERATOR AIR SIDE GAS SIDE 
CMIN/CMA&X CEEPFECTIVENES REYNOLDS NO REYNOLDS A 
o.S67 C.770 4115. 37726 
GS SIDE REGENFRATGR REGENERATCR ATIC GAS/ATI 
DELTA P/P VOLUME WEIGHT RONTAL AREA 
(CE® 13.9924 425.28 1.292 
03 P4 P5 T2 T3 T5 Té 
sieleye 0 | Sir (aa ie Eom | 873.6 1497. 1683. 1079. 
FIGURE 28 COMPUTER OUTPUT 


APPENDIX C 


ERROR ANALYSIS 


It is a basic fact that all engineering investigations are subject 
to imaccuracies. As in any other study, a computer analysis is in- 
herently in error due to: (1) the basic assumptions involved to facili- 
tate the problem solution (as neglecting plate heat transfer resistance 
in equation 4.7) and (2) the approximations made in order to use pre- 
viously compiled engineering data (as the polynomial curve fittings dis- 
cussed in Appendix B). A complete examination of all the errors 
introduced by these factors would be not only time consuming, but quite 
pointless. 

The most significant error in this computer analysis arises from 
the use of tolerance limits on necessary output parameters and finite 
increments on input parameters. This investigation necessitates the 
use of two sets of tolerance levels. The first, Set I, involves a 
Coreétanee™ Wimtt of 3.32°0.3509" = Qe, = 0.5401) om the desired 
thermal efficiency of 0.34 and a limit of 0.2% (599.0 £ HP <£ 601.0) 
on the desired 600 horsepower. Data for these limits is calculated at 
frontal area increments for the recuperator core of 0.02 sq ft and 
length increments of 0.01 ft. Further investigation, in the area of 
optimization, requires tighter tolerances of 1.45% (0.3395 =Qa,4% 0.3405) 
on the thermal efficiency and 0.09% (599.5 = HP = 600.5) on the 
horsepower. Data for these levels, designated as Set II, is supplied 
at frontal area increments of 0.005 sq ft and length increments of 


0.005 ft. Data from Set I is reproduced in Figures 8 through 14 and 
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Tables III through IV; Set II data is presented in Figures 15 through 
22 and Tables V through VIII. 

As there are no performance characteristics available for the 
majority of recuperator cores investigated in this study, no judgement 
can be made as to the errors introduced by the computer simulation. An 
investigation of variations in core parameters as presented in Table VIII, 
however, yields an idea as to the relative error introduced by the neces- 
sary use of finite tolerances and increments. 

It should first be noted in Table VIII that three data sets are 
presented for each core. Set III involves the use of tolerance levels 
as specified for Set II, with reduced frontal area and length increment 
magnitudes of 0.002. Variations between Sets II and III are so small 
that the added computer calculation time introduced with the use of Set 
III over-rides the slight benefits due to the smaller increments. 

Variations in all outputs except dimensional and weight parameters 
are of such an order that the error in computing these quantities seems 
negligible. The most noticeable differences between output sets arise 
in the frontal area and length requirements and the subsequent weight 
and volume quantities. The effect of tolerance levels on dimensional 
requirements is more noticeable in Figures 29 and 30. 

As may be seen, for any given frontal area, the tighter tolerance 
of Set II produces a core of greater length than Set I. A comparison 
of the two plots results in the fact that this effect is more notice- 
able at larger core spacings. For all cores tested, length variations 
between sets decrease as frontal areas increase and less core length is 
required. It is difficult to establish a trend as to the frontal areas 


of minimum volume cores for the two sets. Although Table VIII indicates 


88 


that Set II produces minimum volume cores of greater frontal area than 
Set I, this effect does not hold for all cores investigated. However, 
all Set II minimum volume cores of less frontal area than the corres- 
ponding Set I minimum volume core require a much greater length than 
that for the Set I core. Thus, the overall volume and weight trends 
between cores for each specific set seem to hold. 

The most important questions concerning relative error of the analy- 
sis have been answered. There is evidently little error in the perform- 
ance characteristics output. Dimensional comparisons between the two 
sets result in a maximum difference in length of 0.1 ft for any given 
frontal area and 0.1 sq ft in area for any minimum volume core. Most 
important, due to similar trends in the two sets, the data presented 
may be assumed to provide specific patterns in plate-fin counterflow 
recuperator core dimensions and weight requirements, while the actual 
magnitudes presented may be assumed to be within a reasonable range of 


those which might actually be required. 
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